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Invariant foliations are geometric structures for describing and understanding the qual- 
itative behaviors of nonlinear dynamical systems. For stochastic dynamical systems, 
however, these geometric structures themselves are complicated random sets. Thus it is 
desirable to have some techniques to approximate random invariant foliations. In this 
paper, invariant foliations are approximated for dynamical systems with small noisy 
perturbations, via asymptotic analysis. Namely, random invariant foliations are repre- 
sented as a perturbation of the deterministic invariant foliations, with deviation errors 
estimated. 
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1. Introduction and motivation 

Invariant foliations, as well as invariant manifolds, provide geometric structures 
for understanding the qualitative behaviors of nonlinear dynamical systems, and 
they have been extensively studied for deterministic systems j9l 81 2| . 

Invariant manifolds or foliations for finite dimensional stochastic systems or 
stochastic differential equations (SDEs) were studied in |15l II 12l"5] . Recently the 
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existence of invariant manifolds and invariant foliations for stochastic partial differ- 
ential equations was investigated in |61 71 131 4] and [11] , respectively. In [M] , we esti- 
mated the impact of small noise on invariant manifolds for nonlinear systems. Note 
that random center-like invariant manifolds were approximated for some stochastic 
differential equations (SPDEs) by Wang and Duan [1^, and Blomker and Wang 
[3]. In this paper, we consider a procedure to approximate invariant foliations for 
nonlinear systems perturbed by small noise. We compare invariant foliations for the 
original deterministic dynamical systems and for the randomly perturbed systems. 

We consider the following nonlinear stochastic evolutionary equation with a 
multiplicative noise, in a separable Hilbert space H with a scalar product < •, • > 
and the induced norm 1| • || = \/< •, • >: 

^=AU + F(U) + eUoW, (1.1) 
dt 

where A is a linear (bounded or unbounded) operator, "o" is in the sense of 
Stratonovich stochastic calculus, W = W{t,u!) is a scalar Brownian motion de- 
fined on a probability space (ri,J-',P), and e is a positive parameter representing 
the intensity of the noise. This covers some SDEs and SPDEs. Note that the Ito's 
form of (|l.ll) is 

dU , U , , 

— = A[/ + e- -I- F{U) + e UW. 

The nonlinearity F{U) satisfies F{{)) = and DF{0) = and is Lipschitz continuous 
on H 

\\F{Ui)-F{U2)\\<Lf\\Ui-U2\\, 

where Li? is the positive Lipschitz constant and || • || the norm in the Hilbert space 
H. When the nonlinearity F(U) is locally Lipschitz continuous, the approximation 
result in this paper can be applied to the modified stochastic equation where the 
nonlinearity is appropriately cut-off and thus obtain approximation information for 
the local random invariant foliations. The state space H is the Euclidean space M" 
when the above equation is a SDE or a function space if the above equation is a 
SPDE. When e = 0, Eq. (|l.ip reduces to a deterministic evolutionary equation: 

^=AU + FiU). (1.2) 

We compare the invariant foliations for the original deterministic system (|1.2p and 
for the randomly perturbed system (jl.ip . and quantify their difference when the 
noise intensity e is small. 

This paper is organized as follows. In section 2, we review some basic concepts 
of random dynamical systems, and recall the existence result for random invariant 
foliations. The main result on asymptotic analysis for random invariant foliations 
is described in section 3, and two illustrative examples are presented in section 4. 

2. Random invariant foliation 

Following [7], we assume throughout the paper that the linear operator A : D{A) — > 
H generates a strongly continuous semigroup e"^* on H, which satisfies the pseudo 
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exponential dichotomy with exponents a > > P and a bound K > 0, i.e., there 
exists a continuous projection on H such that 

(i) P"e^* = e^*P"; 

(ii) The restriction e'^*|/j(pu), i > 0, is an isomorphism of the range i?(P") of P" 
onto itself, and we define e"^* for t < as the inverse map; 

(iii) The following estimates hold 

I e^*P"a; | < Ke"* | x |, t < 0, (2.1) 

I e-^^P'x I < isTe^* I a; I, t>0. 

where =I-P''. Denote i/"' = P"iJ, if" = P"if and hence H = H'^H"". 
2.1. Random dynamical systems 

A measurable random dynamical system on Hilbert space (iJ, S) over a driving 
system {9{t))t£T with time T is a mapping 

(fi : T X fl X H ^ H, {t,uj,x) ip{t, uj, x), 

with the following properties [1]: 

(i) Measurability: (p is B{T) ^ (E) ;B-measurable. 

(ii) The mappings Lp{t,uj) = Lp{t,uj, ■) : H ^ H form a cocycle over 0(-), i.e. they 
satisfy (/3(0, cij) = idx for all w G fi and ^p{t + s,uj) = (^(t, d{s)uj) ip{s, uj) for all 
s,t€T and w g 17. 

For SDEs and SPDEs [1], we identify CLi(t) = W{t, uj), and define the driving system 
d{t) is the Wiener shift, i.e., OtUj{T) = lu{t + t) — (^(i). 

To facilitate random dynamical systems study of (jl.ip , we convert it into an evo- 
lutionary equation with random coefficients, called a random evolutionary equation. 
To this end, we introduce z{w) as the stationary solution of the following Langevin 
equation 

dz + zdt = edW. 

Thenz(ii;) = eZ{uj), where Z{uj) is the stationary solution of (iZ(t)+Z(t)(it = dW{t) 
and it can be expressed as 

Z{uj) = / e^'dWir). 

J —oo 

Moreover, 

Jo 

Define a transform 

X ■.= T{uj,X) =Xe-^('^) 

with its inverse transform 

X :=r-i(w,x) =xe^("'). (2.2) 
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Denote U{t,uj,X) as the solution of (|l.ip with initial value X. Introducing 

u = T{etUJ,U{t,Lu,X) = e-'''-^''^'^U{t,Lo,X), 
then the new system state u satisfies the following random evolutionary equation 

m 

du 

— = Au + z{9tUj)u + G(etUJ,u), u(0)=xeH, (2.3) 
dt 

where 
and 

G{uj, u) := e-^('^)i^(e^(")u). (2.4) 

We often denote the solution of (|2.3p to be u = (p{t,uj,x). The solution mapping 
of (|2.3p . i.e. {t,uj,x) — >■ </3(t,a;,a;), generates a random dynamical system. Thus (see 

H) 

{t,Lu,X) T-'^{0tLo,ip{t,io,T{Lu,X)) U{t,Lo,X) 

is also a random dynamical system. In fact, The relationship between solutions of 
(fTTTI) and ([23]) is described by 

Uit, Lu, X) = T-\etio, ifit, cu, T{lu, X)) 
(p{t, oj, x) = T{9tUJ, U{t, oj, T~^{uj, x)). 

2.2. Random invariant foliation 

The concept of invariant foliation is about quantifying certain sets (called leaves 
or fibers) in state space H, starting from all points in such a leaf the dynamical 
orbits have similar asymptotic behaviors. These leaves are thus building blocks for 
understanding dynamics. 

Let us consider a leaf for random invariant foliation for the above random dy- 
namical system (p{t,u},x). A leaf passing through a point $° in the state space 
H, denoted as W{^'^,u!), is a random set and is invariant in the following special 
sense |ll llj 

ip(t,uj,W{<S>°,Lu)) (lW{ip{t,uj,<i>°),etLu) for t>0. 

If we can represent Q as a graph of a C'"' (or Lipschitz) mapping, then (3($°,a;) is 
called a C'^ (or Lipschitz) leaf for the random invariant foliation. The existence of 
random invariant foliation for (j2.3p is shown in |11| . To facilitate our asymptotic 
analysis in the next section, we recall as follows. We only consider stable leaves, still 
denoted as W^($°,w). Unstable leaves may be considered similarly. 
Define 

m = m ~ <&(*), (2.5) 

where l>(t) ~ ip{t,Lo,^^) and $(t) ~ i^(t,a;,$°) are two solutions of (|2.3p starting 
at two initial states $° and respectively. Also introduce the following Banach 
space, for each r], 13 < rj < a, 

C+ = {yj : [0,oo) ^ i7 I ^ is continuous and sup e-'^^^-^o <^-'^)^'^\\ip{t)\\ < oo} 

te[o,oo) 
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with the norm 

M^^= sup e"''*-/o^('^^-)^ni^WII- 

t6[0,oo) 

It is shown in ([H]) that G W{^'^,uj) if and only if there exists a function 
V'(-) e 6+ with V'(0) = and 

Jo 

nt 



where ^ = ^"(4-° - and AG(a;, tA, = G(w, V + *) - G(w, $). Under the gap 
condition 

K^CL f^ + ^) <l, 
there exists an invariant fohation for (12. 3p whose stable leaf is given by 

where <f>° G iJ, $°,Cl)) — > /"(^j $°,a;) is measurable and Lipschitz continuous in 
^ and 

Moreover 

r(^,$°,w) = P"$" + P"V(0;^-P'$",$°,cj), Ce-ff", (2.7) 

where 

rt 



It also follows from ()2.6p that defined in (j2.5p . satisfies the following 

equation 

^=Au + z{9tuj)i^ + AG(0tw, (2.8) 
3. Asymptotic analysis for random invariant foliation 

In this section, we propose an approach to approximate the random invariant foli- 
ation by asymptotic analysis for e sufficiently small. 

Consider the stable leaf of the invariant foliation for (|2.3p (0 < e <C 1), passing 
through a point £ H , 

W{^\u:) = {C + ^\lo)\^ G H^}. (3.1) 

Let the deterministic leaf (i.e. e = 0) be represented as 

U + /('')(0|eei/^}, (3.2) 

where Z"(-, w) : H" H"" and ?('')(•) : ff" are Lipschitz mappings. We expand 

oj) = /('')(0 + el^^\i, $0,6^) + e2/(2)(^^ $0, w) + . . . + e'^7W(C, $°, cj) + • • • 

(3.3) 
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Write the solution of (|2.8I) in the form 

= V'''' {t) + ei/;(i) (t) + • • • + (^^j^^^ {*) + ■■■ (3.4) 
with the initial condition 

m = e + - = e - + /('^)(0 + eZ^')(^,c.) + • • • . (3.5) 

Similarly, write $(<) = (/3(t, w, <f>°) as 

$(i) = $('^)(i) + e$(i)(t) + • • • + e'=$('^(i) + • • • (3.6) 
with initial condition 

$(0) = (3.7) 

By the Taylor expansion, we obtain 

g.(fl.M) ^ ^ezie.i^)) ^ ^ ^ eZiOtii^)) + ■■■+ + . . . (3.8) 

fc! 



and 



/;Z(0,(c^))dr 

1 + ey Z(6l^(Lj))rfr+ ••• + — ^ '— + ■■■. (3.9) 

Suppose F{u) is sufficiently smooth with respect to u. With p.Sp . it follows from 
(1231) that 



G(0tw,?A(t)) = e-"(''*('"»f^ (e"(''*("»i;0(i)) 
= e-^^(«*("»F (e^^('''("))^(t)) 

= (1 - eZiOtiuj)) + . . .) F( (1 + eZ(0t(c.)) + • • •) (v^^^HO + e^^'' W + •••)) 

= F(^('^)(t)) + e (-Z{e,{uj))F{rp'^''\t)) + F,f< U^'^t) + Z{e,{uj))rp'^''\i)]) +■■■ 



(3.10) 



where Fu represents the first order Frcchct derivative [T^ of the function F{u) 
with respect to u and evaluated at 'ip^'''^ (t). In Euclidean space, the Frcchet derivative 
reduces to the classical derivative. 

Substituting p.3p . (|3.10p and p.4p into (|2.3p . and equating the terms with the 
same power of e, we get 



/ = A^W(t) + ^^(V.W(i) + $('^)(t) - F($W(t), 



and 

(it — ^ 

V'(i)(o)-z(i)(e,w), 



^ ^ (3.12) 
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where 



+ pf'>+^^'^ f $(i)(^) + z{9tu)ii;^''Ht) + 



F 



Solve for V^'^HO and V'^HO: 



(3.13) 

= e^V^''HO) + y e-^t*-^) (i^(V'^'^Hi) + ^^'^HO " F {<l>'-'^\t)^ ds, (3.14) 
V,(i)(i) =e^*+^*^'''''^''^^ ('Mi)(e,c.)-y|^-^^+^"^r'^''^'-Arf^^ . (3.15) 



Similarly, we have 



r^£!;M = A$W(t) + F($('^)(i)), 

|$(d)(0) = $0, 



and 



dt ~ 
$(1)(0) 



A + F 



0, 



where 



B = -Z(6lt(cj)) [-$(^)(i) + F($('*)(t))-F*'" $('*)(<) 
Moreover, solve for $('''(f) and $'^^(t), 

$M(t) =e^*4'('*)(0)+ / e'^^*-^)i^($(i)(s))ds. 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



With dlH), dsn and (|3TT0)) . the right hand side of (EH]) can be written as 



P"$° + / e-^'+J'^ "(^'-("^ '^''P"G{e,uj, u(s)) ds = /o + e/i + i?2 
Jo 

where i?2 represents the remainder term and the other two terms are, 

/o=P"<I.o+ /'°e-'*^P" |"f(<I>o(s) + (s)) -F(*o(s))l ds, 
/i = y '^"^"K/ 2(6r(w))dr-Z(04w))j cj ds, 



(3.21) 



with 



(7 =F{i;^''\s) + $('^)(s)) - F(<i>(^))(s) + Ff 

X (V'^^s) + + ^(6'sC^)(V^'''' + $^''^)) . 



(3.22) 
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Substituting p.3|) and p.2ip into (j2.7|) . and matching the powers in e, we get 



I'-'^HO = lo = P''<S>° + 



^-As pu 



F($o(s) + *^''Hs))-F($o(s))l ds, (3.23) 



and 



As a summary, we obtain the following result about approximating invariant fo- 
liation for the random evolutionary equation (|2.3p . including some random ordinary 
or partial differential equations. 

Theorem 3.1 (Approximate invariant foliation for random evolutionary 
equations). Let W((f>°,a;) = + /•'*(^, w, $°) ^ G i/*'} represent a stable leaf, 
passing through a point , of the invariant foliation for the random evolutionary 
equation ^ = Au + z{9tuj)u + G{9tuj, u). Assume that 

(i) F{u) is twice continuously Frechet differentiable with respect to u; 
(ii) For some r/ (a > rj > j3), the following gap condition is satisfied 



K Lf 



< 1. 



(3.24) 



?7 — /? a — rj ^ 

Then for e sufficiently small, the leaf of the random invariant foliation can be ap- 
proximated as 

where \\R2\\ < C(cj)e^ with C{uj) < oo, a.s., 



F(<i>o(s) + vI/W(s))-F(<I>o(s)) 



ds. 



(3.25) 



(3.26) 



4. Examples 

Let us look at two examples. 
Example 1 

Consider a SDE system 

jx = -X + eX oW, 
\y = Y + X^ + eYoW, 

where e > and is a scalar Brownian motion. In this example, A 
H ~ M.'^ {sl finite dimensional Hilbert space), = 



(4.1) 



- 1 
1 

and i?" 
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)w £ M I . The transformed differential equations with random coefficients are 
J 

ix = -x + eZ{etUj)x, 

where Z{ui) is the stationary solution oidZ+Zdt = dW, i.e. Z{uj) = e'^dW{T) 

and ZiOtUj) = e-*Z{io) + e"* e^dWir). 

The stable leaf of the invariant foliation for (|4.2p . passing through a point 
ixo,yo), can be approximated as 

w = {(^)\y-y'> = - ^^^^H - ^^-^ e-^^ diy.) + o{e% . e m} 

Figure [T] compares a random stable leaf (two samples are shown here) for (14. 2p with 
the deterministic stable leaf, passing through the point xo = and yo = 0. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

X 

Fig. 1. The random stable leave passing through the point (0, 0) for Example 1. Two samples of the 
random stable leaf are shown here, together with the stable leaf for the corresponding deterministic 
system (e = 0) 



Example 2 

Consider the following SPDE 

Ut = iU,, + 10U)-U'^ + eUoW,x& [0,1], 
U{0,t) = U{l,t) = 0, 
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where e > and is a scalar Brownian motion. In this example, A = A + 10, 
H = L2(o, 1), D{A) = F|f (0, 1), F{u) = u^. Note that the eigenvalues of A arc A„ = 
10— (nn)'^, and the corresponding normalized eigenfunctions are e„ = sin(n7ra::), 
n = 1, 2 • • • . Here = Span {ei} and H" = Span {e2, 63, • • • , e„, • • • }. 
The transformed random partial differential equation is 

Ut = (mxx + lOu) + Z{6tLo)u - e"2Z(etM)y3^ ^ ^ jg, 1], 
w(0,t) = u(l,i) ^ 0. 

Unlike Example 1, here we can not express Z^^^ (^, w) analytically, but only 
estimate it via p.26p . 
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